Realizing quasi-monochromatic switchable thermal emission from electro-optically induced topological phase transitions

Explorations into the photonic analogs of topological materials have garnered significant research interest due to their application potential. Particularly in planar systems, the prospects of engendering extinguishable topological states can have wide-ranging implications. With an objective of employing these concepts for thermal emission engineering, here, we design and numerically investigate a quasi-monochromatic highly directional mid-infrared source elicited from inversion symmetry-protected topological interface states. Notably, by relying on the architecture of electro-optic effect-induced topological phase transitions, we introduce the possibility of ultrafast switching of thermal radiation. These reversible phase transitions, being free from carrier transport are inherently fast and evoke thermal emission modulation with a modulation depth upto 0.99. Specifically, our platform exhibits a near-perfect extinguishable spectral emission peak at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4~\mu$$\end{document}4μm with a quality factor of over 18500, displaying negligible parasitic emissions. Furthermore, the optimized interface state manifests itself for only one of the polarization modes, resulting in polarized emission under resonance conditions. To establish a methodical approach to parameter optimization, we also model our platform as a leaky mode resonator using the framework of temporal coupled-mode theory. We believe, our findings can provide a way forward in establishing complete control over the optical characteristics of the infrared thermal emitters.

Realizing quasi-monochromatic switchable thermal emission from electro-optically induced topological phase transitions Nitish Kumar Gupta 1* , Sapireddy Srinivasu 1 , Anjani Kumar Tiwari 2 , Harshawardhan Wanare 1,3 & S. Anantha Ramakrishna 3,4 Explorations into the photonic analogs of topological materials have garnered significant research interest due to their application potential. Particularly in planar systems, the prospects of engendering extinguishable topological states can have wide-ranging implications. With an objective of employing these concepts for thermal emission engineering, here, we design and numerically investigate a quasimonochromatic highly directional mid-infrared source elicited from inversion symmetry-protected topological interface states. Notably, by relying on the architecture of electro-optic effect-induced topological phase transitions, we introduce the possibility of ultrafast switching of thermal radiation. These reversible phase transitions, being free from carrier transport are inherently fast and evoke thermal emission modulation with a modulation depth upto 0.99. Specifically, our platform exhibits a near-perfect extinguishable spectral emission peak at 4 µ m with a quality factor of over 18500, displaying negligible parasitic emissions. Furthermore, the optimized interface state manifests itself for only one of the polarization modes, resulting in polarized emission under resonance conditions. To establish a methodical approach to parameter optimization, we also model our platform as a leaky mode resonator using the framework of temporal coupled-mode theory. We believe, our findings can provide a way forward in establishing complete control over the optical characteristics of the infrared thermal emitters.
Attributable to its statistical nature, the thermal emission is considered to be lacking the features of selectivity and hence almost exclusively associated with characteristics such as broadband spectrum, unpolarized emission, and a quasi-isotropic angular emittance 1 . Therefore, the prospects of exhibiting control over thermal radiation are of widespread fundamental significance, apart from being crucial to the design and operation of devices ranging from micron-scale chips to astronomical size systems. Inter alia, the fruition of thermal emission with ultranarrow bandwidths has recently garnered significant research interest owing to its direct utility in applications ranging from -high figure-of-merit infrared (IR) sensing, IR spectroscopy and fingerprinting, thermophotovoltaics, photon mediated thermal logic, IR labeling, anti-counterfeiting, and in countermeasures development against IR homing projectiles 2-6 . The development of techno-commercially plausible customized thermal sources, particularly in the mid-IR domain, however, has its challenges. Owing to the fact that the spontaneous emission scales at cubic power of frequency, an IR analog of LED turns out to be highly inefficient 7 . Although state-of-the-art solutions like Quantum Cascade Lasers can surpass most of the spectral purity requirements, the intersubband transitions based operation requires stacking up of many hundreds of quantum wells and barrier layers 8 . Motivated by these viability concerns and to develop application-specific minimalistic solutions, many researchers have revisited the problem of optimizing the radiant emittance per unit area from an object (defined in terms of the Planck's law of spectral radiance) 9 and noticed that although the temperature management is one obvious way to control the thermal emission; however, the theory also suggests a trickier but more rewarding alternative of emissivity engineering [10][11][12][13][14][15][16][17][18][19][20] . This realization has brought the nano, and microphotonic structures to the forefront, which can artificially alter the absorptivity/emissivity using structure-derived electromagnetic www.nature.com/scientificreports/ resonances. Such designs provide an austere yet amply suited template to control the light-matter interaction in a way that near-perfect spectrally selective, efficient, and highly directional absorption can be materialized. As established by Kirchhoff 's law, under thermodynamic equilibrium, the reciprocity can be factored in and the same concepts also provide an equivalent assessment of emission. Over the years, however, it has been observed that many of these proposals also resort to complicated lithography techniques 1, 21 , thereby raising concerns of structural stability at high temperatures apart from a limited commercial acceptance. Additionally, significant background and multiple resonances lead to high noise levels and wasteful parasitic emissions, affecting the conversion efficiencies and the achievable signal-to-noise ratios. Furthermore, these proposals cater only to the static spectral selectivity requirements. High-speed switching and modulation have been perennially considered incompatible with thermal emission due to inherently slow dynamics of heat transport, even for devices with small thermal inertia. This had limited the modulation speeds up to 100 Hz 7,22 ; only of-late the proposals have been put forward to achieve MHz switching [23][24][25] . However, a thermal emission modulation scheme free from carrier transport and complying with thermodynamic equilibrium can offer orders of magnitude faster speeds, motivating us to look for atypical replacements. On a different note, the application of topology, a concept of modern mathematics, to the condensed matter physics led to the discovery of neoteric and exciting paradigms such as topological phases of matter and effective magnetic fields 26 . Particularly, in Fermionic systems, the inception of band topology dates back to the seminal work by Thouless et al. 27 , which has been gradually extended to classical electromagnetic wave systems 28 . The keystone of the formalism lies in conception that the momentum space global properties of a topologically nontrivial bulk bandstructure can be encapsulated in quantized topological invariants 26 , and by combining structures of different topological orders, unconventional interface or edge states can be manifested with detailed control. These topological interface states (TIS) provide a doorway to additional degrees of design freedom, consequently garnered enormous research interest and facilitated the dawn of topological photonics 26,[28][29][30][31][32][33][34][35] .
Our present work concentrates on employing the principles of topological photonics to conceptualize topological phase transitions in a spatial inversion symmetric (SIS) multilayer configuration, which being relieved from inertial aspects, can lead to the ultrafast switching of TIS. The dynamic topological phase transition in our structure is induced by the repetitive system evolution across a Brillouin Zone (BZ) edge bandgap closing point using the electro-optic (EO) effect of anisotropic LiNbO 3 films. The resulting symmetry-protected TIS is then integrated with a localized loss of atomically thin monolayer graphene, leading to a monochromatic thermal emitter with temporally agile modulation capabilities. Notably, the movement across the topological phase transition point assigns a new degree of design freedom to the realized interface state where its very existence can be controlled by weak modulations in the refractive index. The platform is designed to operate in a technologically important atmospheric high transmission window at mid-IR ( 3 − 5 µm). Specifically, the optimizations have been performed to synthesize the emitter at the middle of a photonic bandgap at 4 µ m; this provides high selectivity and rules out the possibility of parasitic emissions leading to high conversion efficiencies. Furthermore, by circumventing the patterning requirements, our platform not only ensures structural stability at high operating temperatures but also provides means for mass production using standard deposition techniques.

A routine for bandstructure engineering
The platform under consideration here consists of a binary one-dimensional photonic crystal (1D PhC) constituting from spatial inversion symmetric unit cells. In such a stratified medium, the Bloch waves E z n,K (x) = u n,K (x)e −iKx define the eigenmodes of the system, with u n,K (x) being the periodic part of the Bloch solution, matching the periodicity of the structure , i.e. u n,K (x + �) = u n,K (x) . The bandstructure calculations for this system can be performed using the dispersion relation 36 where K is the Bloch momentum; Z i = (µ i / ǫ i ) is the impedance and k i is the wavenumber in the i th medium ( i = a, b ); is the lattice constant ( = d a +d b ). The absolute value of the right-hand side (RHS) of the Eq (1) assigns the frequency regions corresponding to the bands and bandgaps. When it is smaller than or equal to 1, the equation gets satisfied for a real value of Bloch momentum K, and the corresponding solutions are the propagating modes defining a passband. On the contrary, if the absolute value of the RHS is greater than 1, only imaginary values of K can satisfy the equation leading to exponentially decaying states denoting a bandgap. It has been explained earlier 37 that the maxima and minima of RHS define the midgap frequencies, whose location can be approximately ascertained from the expression ω m = m πc n a d a +n b d b (with m being an integer and n a , n b are the refractive indices corresponding to the medium 'a' and 'b'). Therefore, although the bandstructure itself is sensitive to the unit cell composition, the positions of midgap frequencies will remain oblivious to it, determined solely by the total optical thickness ( n a d a + n b d b ). This provides an opportunity for calibrated manipulation of the bandstructure and a venue to study bandgap/band transitions. To fathom the resulting implications, we proceed in a systematic manner by first quantifying the fractional composition of the unit-cell (depicted in Fig. 1a) and then using it to study the parametric performance of a representative PhC. Accordingly, we define the fractional weightage of material 'a' in the unit cell as α = n a d a n a d a +n b d b and parametrize the individual layer thicknesses in terms of it. Thereafter keeping the optical thickness constant, we scan the entire parametric space of α to obtain the reflectance performance of a PhC consisting of 10 periods with n a = 2.0 and n b = 3.0 (the design wavelength corresponding to 1 st order bandgap is taken to be 4 µm ). These calculations are performed using the standard transfer matrix method (TMM), and the results are plotted in Fig. 1b, highlighting the following attributes: (1) from the topological standpoint, the first order bandgap is trivial as it entails no prospects www.nature.com/scientificreports/ for band inversion; only from 2 nd bandgap onwards we see the occurrence of probable band inversion sites. (2) As we move to higher-order bandgaps, the fractional change required in the unit-cell composition for an openclose-open transition comes down significantly. Moving on, it is noteworthy to mention that if we are interested in realizing a topological interface state only, this route of thickness parametrization turns out to be more conducive for devising a fabricable design and will be implemented elsewhere. However, in this work, since we aim for temporal modulation of TIS, it warrants dynamic manipulation of topological invariants of the system. Therefore, external means of dynamic bandstructure tuning must be called in. To this objective, we resort to the EO effect-induced perturbative changes in the refractive index (RI). For exploring the sites of band inversion, we look for refractive index parametrization in a representative PhC. Such a parametrization, however, is not as forthright and requires more careful investigation as a change in the RI contrast would have wide-ranging impacts on bandstructure. Nevertheless, we begin by defining a new parameter ρ = δn i /n i which quantifies the fractional change in RIs of PhC constituents around their base values (keeping optical thickness constant) and obtain the reflectance characteristics of a 10 period PhC as done previously (with base RIs n a = 2 , n b = 3 and fixed thicknesses of d a = 0.5µm , d b = 0.33µm corresponding to α = 0.5 ). We also want to point out that the configuration adopted here would lead to an increase in RI contrast with an increasing value of ρ (which is more aligned with the case of LiNbO3-Si platform to be used later on). Notice here that the parametric space of ρ is highly restricted not only on account of the practical limitations but also keeping in mind that it would make little sense to talk of ρ = 1 . The calculation results are plotted in Fig. 1c, sufficiently highlighting the fact that the small RI alterations possible with the EO effect will always be insufficient in realizing the band inversion for lower-order bandgaps.
The ongoing analysis, therefore, encourages us to safely assume that in the present work, we have to invariably work with higher-order bandgaps. Equipped with these insights, in the next section, we realize a TIS modulator around the 5 th order bandgap for a concatenated PhC structure.

Results: electro-optically switchable topological phase transitions
Our objective is to first realize a PhC platform admitting an EO effect-induced manipulation in its bandstructure (and a possible band inversion). We commence the design with material considerations: on account of the proven compatibility of LiNbO 3 with the silicon-based platforms, the constituent materials have been selected to be Si and LiNbO 3 having a RI of 3.45 and 2.04 respectively in mid-IR 38,39 . The electro-optically induced change in LiNbO 3 RI can be ascertained by the following relation: where �( 1 n 2 ) ij is the second-rank tensor of the change in relative permittivity, r ijk is a third-rank linear EO coefficient tensor, E k is electric field component along k; i,j,k correspond to x,y, and z coordinate axes respectively. Among all the EO coefficients r ijk , the r 33 is the largest linear EO (Pockel's) coefficient (see Supplementary  (1), within the mentioned span of available RIs, for three representative scenarios have been plotted in Fig. 2a: Case-1 marks the reference case under no external field with Si and LiNbO 3 thicknesses of 3.08 and 1.07 µ m; Case-2 corresponds to a moderate application of field and the associated reduction in LiNbO 3 RI to 2.00; while Case-3 denotes a further reduction to a value of 1.94. Noteworthy to mention that in order to illustrate the concepts, thicknesses in Case-2 and 3 have been adjusted for ensuring the invariance of band centers. Thereafter, the topological character of the bulk bands has been ascertained by the extraction of geometric Zak phases 37 calculated using Eq (3) and denoted alongside the bands in bold black letters.
where θ zak n is the Zak phase associated with n th band; the quantity in the parenthesis is the Berry connection (or Zak connection) of the n th band ( A Zak n ); ǫ(z) is spatial permittivity function and u n,K (z) is the Bloch eigenfunction in the n th band, which is obtained using the established eigenvalue formulation followed in stratified media 36 . From the bandstructure evolution presented in Fig. 2a and b and the corresponding Zak phase calculations, we notice that the open-close-reopen transition at the 5 th order bandgap indeed corresponds to a band inversion and an associated reversal of topological invariants (Zak phases) for the neighboring bands. This proclaims a transition in momentum space topology of bulk bands with a change in LiNbO 3 RI . The evolution of bandgap width and Zak phases have been specifically captured in Fig. 2c and d with RI contrast (defined as difference between the RIs of Si and LiNbO3). According to the bulk-interface correspondence, ramifications of such a transition should be traceable as a sign reversal in the surface impedance of the associated bandgaps 37 , which can be ascertained by following the evolution of Zak phases of all the bulk bands below it, resulting in www.nature.com/scientificreports/ where ξ (n) is defined as iξ (n) = Z (n) Z 0 ; l is the number of crossings under the n th bandgap. Calculations as per Eq (4) confirm this proposition which has been explicitly represented in Fig. 2a by color-coding the bandgaps. In this formalism, we can see that the surface impedance of the 5 th order bandgap reverses its sign (color) during the transition from Case-1 to 3, engendering the possibility of a topological interface state.
To witness the manifestation of such a TIS, we create a defect-free assembly of PhC-I and PhC-II (in the sense that there is no defect layer incorporation and no break in periodicity) having periods N 1 = N 2 = 10 , and obtain the concatenated structure of Fig. 3a, which will be termed as a TIS modulator. While the PhC-I always retains the character of Case-1, the PhC-II can make a phase transition from Case-1 to Case-3, depending upon the externally applied field and the resulting perturbative change in RI. During this transition, the surface impedance sign reversal translates to reflection phase reversal as plotted in Fig. 3b. It brings us to a dynamic conjugate impedance matching condition: when the modulator is in 'ON' state, a TIS emerges at 4 µ m (Fig. 3c) on account of mismatch in topological invariants of the two constituent PhCs. In the 'OFF' state, however, the TIS gets extinguished (see Fig. 1 in Supplementary Information) as the topological invariants of PhC-I and PhC-II become equal. For vindicating the interface state nature of TIS, the electric field magnitude profile of TIS at the resonant wavelength of 4 µm has also been provided in Supplementary Information.
Before moving any further, we also succinctly explain the polarization-related modalities of the obtained TIS. In the concatenated PhC structure of Fig. 3a, we have chosen the optic axis and the applied voltage to be along the z-axis (with light propagating along the x-axis). As noted before, in such a configuration, the Pockel's coefficient r 33 comes into the picture for z-polarized light, and the corresponding maximum change of the RI in n z is given as where n e is extraordinary RI. The applied electric field instigates a RI change in the n y as well, which can be determined as δn y = − n 3 o 2 r 13 E z where n o = 2.11 is ordinary RI at = 4 µm . The EO coefficient r 13 employed here is about five times smaller in value than r 33 in the wavelength region of interest; in consequence, the RI change δn y would be sufficiently smaller than δn z . This disparity will ensure the occurrence of topological phase transition for only one of the polarization modes (z-polarization) and bring polarization discrimination into the structure. www.nature.com/scientificreports/ With the architecture of a switchable TIS in place, we introduce localized ohmic loss to the structure in the form of a monolayer Graphene 43 (Fig. 4a), the optical properties of which have been modeled using the Kubo formalism 44 , accounting for both, the interband and intraband contributions to the graphene conductivity (see Supplementary Information). The arrangement is then optimized for N 1 and N 2 to realize near-perfect absorption/emission with results plotted in Fig. 4b for N 1 = 15 and N 2 = 16 . Manifestation of TIS in ON-state evinces out an ultra-narrow spectral resonance with an emissivity peak of 0.97 at 4 µ m, leading to coherent infrared emission. The corresponding linewidth and associated quality factor come out to be 0.0002 µ m and 18600. The OFF-state response is also plotted in red, demonstrating a high modulation depth (0.99) and extinction ratio ( −19.4 dB) for our platform. The associated scattering performance is plotted in Fig. 4c with black curve denoting a near-perfect reflection minimum. The electric field magnitude profile for this optimized thermal emitter has also been plotted in Fig. 4d at the peak emission wavelength, which unequivocally illustrates the enhanced energy density at the location of localized loss element (monolayer graphene). Furthermore, as alluded to previously, the conditions of existence of TIS will be satisfied for only one of the polarization modes; thereby, the thermal emission emerging from the structure would be inherently polarized (see Fig. 3 in Supplementary Information).
Although the optimized design of Fig. 4a can attain near-unity emissivity values, analytical curiosities to obtain perfect emission can also be facilitated. To systematically accomplish this and to gain crucial insights into the operating mechanism, we resort to the theoretical framework of temporal coupled-mode theory (TCMT) 45 , which captures the behavior of our thermal emitter as a leaky mode optical resonator. As per the TCMT framework, a single-mode optical resonator coupled to the outside world through m-ports can be modeled with a symmetric Lorentzian profile with absorption around the resonance frequency ( ω ), given as where γ rad is far-field radiative damping rate; γ abs is intrinsic loss rate due to ohmic damping mechanisms; m is the number of ports.
Whilst the device architecture in our platform suggests the number of ports (m) to be two, for attaining perfect absorption we reduce it to single port by blocking transmission . This is accomplished by realizing near-perfect Bragg-reflection from the PhC-II (by setting N 2 = 30 ), which also prohibits the transmission across the platform. Under such conditions, the absorption in the resonator simplifies to 4γ rad .γ abs (γ rad +γ abs ) 2 . Evidently, the key to achieving perfect absorption (efficiency maximization) is to establish a balance between γ rad and γ abs , leading to the critical www.nature.com/scientificreports/ coupling (CC) of incoming radiation. It is noteworthy that the condition of critical coupling stands independent of the strength of the intrinsic loss channel. Hence, even vanishingly small amounts of loss can lead to perfect absorption. This observation provides an alternative analytical explanation to prior results of near unity emissivity in our platform (for a freestanding monolayer, Graphene absorption is only 0.023). The absorption bandwidth, on the other hand, is dictated by the strength of all damping mechanisms present in the system (γ rad + γ abs ) . The many implications and predictions of the TCMT framework in a reduced dimension 2D abstract parameter space ( γ rad , γ abs ) have been summed up in Fig. 5a. The order of the associated damping rates has been ascertained by fixing the operating point to the working wavelength of 4 µ m. While the efficiency of absorption resonance is dictated by the proximity of an operating point to the CC contour, the resonance bandwidth (and the Q-factor) exhibits a monotonically decreasing (increasing) behavior along the marked direction. The contiguous parameter space is then searched for a local maximum with the corresponding running variable in the device to be the number of periods (N 1 ) in PhC-I. The mapping of variation in N 1 to γ rad and γ abs has been done by fitting the TCMT responses to the TMM calculations, and the results are presented in Fig. 5b. Starting from N 1 = 12 , as we increase N 1 , the operating point moves from overcoupled to undercoupled regime, crossing the critical coupling contour at N 1 = 16 , which marks the TCMT prediction for perfect absorption. An exhaustive parameter search of Fig. 5c performed solely based on TMM analysis provides testimony to such claims. From spectral characteristics, we now move on to the spatial characteristics of the designed thermal emitter. It is noteworthy to ponder upon the directional emission characteristics of the proposed thermal emitter, which is ascertained from the angle-resolved emissivity pattern at the operating wavelength, as depicted in Fig. 6a. From these patterns, the angular width of the radiation lobe comes out to be 27 mrad, establishing its ultra-narrow angular emission behavior. In Fig. 6b, we plot the angular width of the emission as a function of N 1 . It depicts that the FWHM of the emission lobe exhibits a monotonically decreasing behavior with an increase in the number of periods, and hence, the directional emission characteristics can further be improved by incorporating more periods in the platform. The value of peak emission, though, does not exhibit such a monotonic behavior and, in line with the predictions of the TCMT framework, starts to fall off beyond an optimum point.
Lastly, we mention some design considerations from the viewpoint of fabrication conduciveness. Firstly, the structure in Fig. 4a can be simplified to some extent by replacing the monolayer graphene sheet with an absorbing film. Such a design would entail a compromise in thermal emission performance owing to the relatively distributed nature of the loss element, but can be instrumental in situations demanding fabrication simplifications. Although the non-negligible thickness and the resulting modifications brought about by a thin absorbing film do not hamper the existence of TIS; however, as established by the TCMT framework, the critical coupling www.nature.com/scientificreports/ condition would be altered, leading to degradation in spectral emission characteristics. Specifics and the calculated performance for one such design incorporating thin plasmonic metal film (of gold) are provided in Fig. 4 in Supplementary Information. Therefore, although graphene is used only as a loss element in our structure and has no role to play regarding the qualitative response of TIS, its incorporation becomes critical for obtaining high quality factors. An even greater fabrication simplification is possible by resorting to a corresponding in-plane configuration, as provided in 46 . The nanobeam PhC configuration discussed in this work deviates from the vertical stacking of layers; instead, the constituent equivalent layers are stacked in-plane with a corresponding change in the light propagation direction.
In conclusion, we have designed and analyzed a quasi-monochromatic and highly directional thermal emitter by harnessing a topological interface state. The TIS owes its existence to the EO-induced small perturbations in the LiNbO 3 RI and the associated modalities of topological phase transitions, which being free from carrier transport mechanisms and heat diffusion dynamics can facilitate ultrahigh-speed mid-IR signal processing. Our approach of capitalizing on a defect-less and lithography-free geometric phase-based architecture offers advantages of volume production and ensures structural stability, opening the prospects for the unobtrusive development of mid-IR integrated photonics hot thermal emitters. www.nature.com/scientificreports/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.